NMR quantum computation with indirectly coupled gates 

David Collins/ K. W. Kim/'0 W. C. Holton/ H. Sierzputowska-Gracz,^ and E. O. Stejskal^ 

^ Department of Electrical and Computer Engineering, Box 7911, North Carolina State 

University, Raleigh, North Carolina 27695-7911 

"^Department of Biochemistry, Box 7622, North Carolina State University, Raleigh, North 

Carolina 27695-7622 
^Department of Chemistry, Box 8204, North Carolina State University, Raleigh, North Carolina 
Q> ; 27695-8204 

gN ; (February 1, 2008) 

-(— > 

o 

Q ; Abstract 

An NMR realization of a two-qubit quantum gate which processes quantum 
^ ■ information indirectly via couplings to a spectator qubit is presented in the 

^^ , context of the Deutsch- Jozsa algorithm. This enables a successful comprehen- 

^D I sive NMR implementation of the Deutsch-Jozsa algorithm for functions with 

three argument bits and demonstrates a technique essential for multi-qubit 
On ' quantum computation. 
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Nuclear magnetic resonance (NMR) spectroscopy has emerged at the forefront of experi- 
mental quantum computation investigations Il|-10|. Key concepts such as fundamental gates 



|]1|J^Q] and error correction [^ have been demonstrated using NMR spectroscopy. However, 
comprehensive algorithm realizations have only been accomplished for the Deutsch-Jozsa al- 
gorithm for functions with one and two bit arguments |lf-p[ and for Grover's algorithm with 



a two bit register ||lTHr3| . In these instances crucial two-qubit gates were realized "directly" 
via the coupling between the corresponding nuclear spins. For a quantum computer with a 
larger numbers of qubits, the associated requirement of appreciable coupling between any 
pair of spins raises difficulties. First, it may not be possible to find any molecule with this 
coupling configuration. Second, it demands increasingly complex schemes for managing the 



evolution of spectator spins during execution of any two-qubit gate ||lj] . However, efficient 
"indirect" realization of any two-qubit gate via a chain of couplings through intermediate 
spins is possible [|15|. This relaxes the coupling requirements and a "linear" coupling con- 
figuration, for which the pattern of couplings is A — B — C — D — ... suffices for quantum 
computation. In this Letter we report a comprehensive three qubit NMR realization of the 
Deutsch-Jozsa algorithm for functions with three bit arguments, using indirect realizations of 
two-qubit gates and a linear coupling configuration for information processing. The method 
that we present is general and readily scalable to larger numbers of qubits. 

The Deutsch problem |lT6|jr7| considers / : {0,1}^ — > {0,1} that are constant or bal- 
anced. A balanced function returns as many times as 1 after evaluation over its entire 
range. Given any / which is either balanced or constant the problem is to determine its 
type. For classical algorithms that solve the Deutsch problem with certainty the number of 
evaluations of / grows exponentially with A^. A quantum algorithm requires a single evalu- 
ation of / and yet solves the problem with certainty [1^,17]. Previous NMR demonstrations 



|I|-0 used the Gleve version [|T^], requiring an A^ qubit control register for storing function 
arguments, plus a 1 qubit function register for function evaluation. Our recent modification 
to the algorithm (see Fig. |I|) only needs the A^ qubit control register [|TB|. It was also shown 



that for each admissible (i. e. constant or balanced) function for N < 2 the evolution step is 
a product of single qubit operations. For an isolated spin | nucleus any single qubit opera- 
tion amounts to a rotation of the magnetization vector; this can be implemented classically. 
Thus for N < 2 the algorithm can be executed classically. Previous comprehensive NMR 
implementations of the algorithm fall within this classical regime []1]-§1 . However, for A^ > 3 
there exist balanced functions for which qubits in the control register become entangled; 
this is indisputably quantum mechanical. A^ = 3 is then the critical point at which quantum 
mechanical features become essential in the Deutsch-Jozsa algorithm. 

Two-qubit, entangling operations only appear during the function evaluation step, Uf. 
However, for / constant. Iff = I and the algorithm can be executed classically. To assess 
quantum behaviour the corresponding gate for balanced / must be investigated. Balanced 
functions can be characterized by the choices of 2^~^ arguments out of 2^ possibilities for 
which is returned. For A^ = 3 this gives f ^) = 70 balanced functions. Admissible functions 
may be represented via power series expansions in the argument bits, Xj where i = 0,1,2. 
That xf = Xi for Xi G {0, 1} implies that for any admissible function, 

2 2 

f{x2, Xi, Xo) = Yl ^ij^i^j ® H ^i^i ® " (l) 

i>j>0 i=0 



where addition is modulo 2 and aij,ai,a G {0,1}. The disappearance of a cubic term in 
Eq. (|lD is a property of balanced and constant functions. This provides a preliminary 
decomposition: 

Uf = li {u'T n {uT (2) 

j>i>o fc=o 

where 

U'^ \x) := (-l)^'^^^ \x) , and (3) 

U^ \x) := (-1)"' |a;) . (4) 

are the quadratic term gate and linear term gate, respectively. The constant term merely 
provides an identity operation. Quadratic and linear term gates all commute and can be re- 
arranged at will. In terms of fundamental gates, W = -R^(180) where superscripts index the 
qubit to which the 180° single qubit rotation is applied. Similarly, W^ = R-'_^{90)U^j^R-'^{90) 

where U^^^ is a controlled-NOT gate with control i and target j. Thus the algorithm can 
be implemented classically for functions with no quadratic terms. However, a quadratic 
term gate can produce entangled states from unentangled states for qubits i and j via its 
constituent controlled-NOT gate. There are no classical operations involving magnetiza- 
tion vectors of two spins that produce entanglement. Therefore, for A^ = 3 it is in the 
quadratic term gates that quantum mechanical features appear. The arrangements of linear 
and quadratic term gates that constitute Uf for any admissible / can be classified via simi- 
larity under permutations of control register qubits. Equations [(0)-(^] imply that this can 
be accomplished by classifying admissible functions via similarity under argument bit per- 
mutations and/or addition of a constant term. Accordingly there are ten classes of balanced 
functions; a representative of each is provided in Table |I[ Function evaluation steps for any 
members of a given class differ only by a permutation of the control register qubits to which 
they are applied. Thus a realization of the A^ = 3 Deutsch algorithm is comprehensive if the 
algorithm is applied to at least one representative from each class of admissible functions. 
All possible quadratic term gates are required for the classes represented by /g and /iq. 
Therefore, for A^ = 3 the algorithm requires two- qubit gates between all pairs of qubits and 
is suitable for testing quantum information processing with a linear coupling configuration. 
Currently the most accessible experimental technology for quantum computing is NMR 
spectroscopy of spin ^ nuclei of appropriate molecules in solution [|T]-0. Any molecule con- 
taining three distinguishable, coupled spin | nuclei in an external magnetic field provides 
the three qubits needed to solve the A^ = 3 Deutsch problem. To a good approxima- 
tion the Hamiltonian for a room temperature solution state sample is i^ = I]^=o ^^l + 
■^J2'i>j>oJijO'ld'i, where cUj are the Zeeman frequencies, Jij the scalar coupling constants, 
and (J* Pauli operators |T9|. Superscripts label the spins and identify them with the corre- 



sponding argument bits. A literal translation of the algorithm into NMR operations would 
begin with initialization via a pseudo-pure state preparation scheme 0-§[. The evolution 
stage can be implemented by building single qubit and controlled-NOT gates from standard 
sequences of spin selective RF pulses and periods of evolution under H [^J^. Measurement 
can effectively be accomplished via tomography, which requires repeated execution of the 
initialization and evolution stages, each time varying the readout pulse before acquisition [HI . 



It is, however, possible to apply the evolution stage directly to a thermal equilibrium initial 
state and successfully solve the Deutsch problem with an expectation value measurement 
PJ20[] . Analysis of the system state after the function evaluation step demonstrates this for 



the A^ = 3 Deutsch algorithm. Using the product operator formalism, the deviation part 
of the thermal equilibrium equilibrium density operator for a weakly coupled homonuclear 
NMR system is proportional to pth = -^^ + -^i + -^° [|H]- The initial rotation with phase 
= — I transforms this to —/J — I]. — /°. Then tff produces the states listed in Table |. 
Signal acquisition immediately after application of Iff and with no additional readout pulses 
provides a spectrum, the f -spectrum, whose form depends on /. A fiducial spectrum is ob- 
tained in the same fashion but with f/j replaced by J; here the system's pre-acquisition state 
is — JJ — I^ — Ix- Comparison with pf for admissible functions (see Table |ip indicates that for 
each there is either a or a vr phase difference between each line each line of the /-spectrum 
and its counterpart in the fiducial spectrum. More precisely for applicable functions: 

(i) / is constant if and only if the f -spectrum is identical to the fiducial spectrum and 

(ii) / is balanced if and only if there is a ir phase difference between at least one line of 
the f -spectrum and its counterpart in the fiducial spectrum. 

This criterion requires that each spin is coupled to at least one other spin. If any spin 
is completely uncoupled then the entire /-spectra for /g and /lo disappears; the comparison 
would be impossible. However, if each spin is coupled to at least one other then for A^ = 3 at 
least one is coupled to the other two. This ensures that at least one of the doubly antiphase 
multiplets for /g and /lo survives, giving a line in the /-spectrum with a tt phase difference 
relative to its fiducial spectrum counterpart. For /y and /s at least one of the antiphase 
multiplets (for spin 2 or spin 0) must survive; again this provides a line whose phase differs 
by TT. For /4, /s and /e the entire spin multiplet in the /-spectrum displays a vr phase 
difference relative to its fiducial spectrum counterpart. The same is true for spin 2 in the 
cases /i, /2 and /s. 

The fiducial spectrum can be phased so that its constituent lines all appear upright. 
Thus, the answer to the A^ = 3 Deutsch problem may be determined by inspecting the 
/-spectrum for inverted lines. Each balanced function produces at least one inversion. 
For constant functions all lines are upright. This provides a solution state NMR scheme 
for conclusively answering the A^ = 3 Deutsch problem with just one application of the 
evolution stage (here equivalent to the unmodified version followed by a i?n(90) readout) to 
the thermal equilibrium input state. 

A saturated solution of ^'^C labeled alanine in D2O provided the qubits. We label the 
carboxyl carbon, spin 2, the a carbon, 1 and the methyl carbon, 0. Protons were decoupled 
using a standard heteronuclear decoupling technique. Scalar couplings are J21 = 56Hz, 
Jio = 36IIz and J20 = 1.3IIz. Relaxation times are Ti(2) = 11.5s, T'i(l) = 1.2s, and 
Ti(0) = 0.7s and T2(2) = 1.3s, T2(l) = 0.41s, and r2(0) = 0.81s where the argument labels 
the spin. Spin selective rotations were implemented via Gaussian shaped pulses of duration 
0.7ms for spins and 1 and 0.5ms for spin 2. No hard pulses were used. Linear term gates 
can be implemented via spin selective phase shifts on the output spectrum by placing them 
after the quadratic term gates and immediately prior to acquisition. Thus certain /-spectra 
differ by spin selective phase shifts only. These are: (i) /i, /2 and f^, (ii) /4, /s, and /g, 
(iii) /7 and /§ and (iv) /g and /lo- The crux of the experiment is in the realization of 
the quadratic term gates. This is accomplished via the rotation and delay construction of 



controUed-NOT gates 0,0,||. Selective refocusing sequences [pl|-p3| effectively eliminate all 
but one coupling term in H, thus providing appropriate evolution during the delay. The 
resulting quadratic term gate simplifies to 

W^ = [1/2J,,]^^- - m\ - [90]i, (5) 

where [9]'^ indicates a rotation of spin i about the axis n through angle 9 and [t]*"' evolution 
under the scalar coupling between spins i and j for period t. For alanine this is satisfactory for 
f/^^ and f/^°. However, for [/^° it is inadequate since 1/2 J20 = 0.42s which is comparable to 
the smallest T2. An alternative is to process the information via spin 1 and use only the linear 
coupling configuration ( spin 2 - spin 1 -spin 0). An indirect realization is f/^° = ^sw^cn^sw' 
where t/°^ is the SWAP gate between qubits 1 and 0. After simplification, 

tj^' = ml - ml - [i/2j,or - mi - ml 
-[i/2j,,r-mi-[i/^j2ir-mi 
-[i/2j,,r-mi-mi-[i/2j,,r 
- ml - m-y - [90]!. - [90]-!, . (6) 

This gives a pulse sequence of duration 0.071s (excluding z rotations that are equivalent to 
phase shifts in the output spectrum) that is significantly faster than that using [I/2J20] . 
During [1/2 Jio] and [1/2 J21] evolution periods the selective refocusing scheme effectively 
removes the scalar coupling between spins 2 and 0. Throughout this implementation of 
the algorithm information is processed using only the spin2 - spin 1 -spin linear coupling 
configuration and not the spin 2 -spin coupling. However, the latter must be taken into 
account in the interpretation of the output spectra. 

The experiments were performed at room temperature using a BRUKER 500-DRX spec- 
trometer and an inverse detection probe. For each representative function listed in Table 
I signal acquisition takes place immediately after implementation of JJf. Figure ^ provides 
selected experimental spectra that are phased so that — /^. product operator terms corre- 
spond to upright multiplets. The line orientations agree with those predicted from pj and 
provide correct solutions to the Deutsch problem. An estimate of errors for the most com- 
plicated case, /g, may be obtained by applying a selective 90° readout pulse about the x 
axis immediately after Uj. Ideally the readout spin multiplet should remain while the others 
disappear. The average amplitudes of the residual signals for the latter lie between 14% and 
31% of the average amplitude of the corresponding lines with no readout. The ability to 
extract the Deutsch problem solution for A^ > 3 via pure phase information in the output 
spectrum, in contrast to amplitude information, appears to have mitigated such errors. It is 
not yet clear how this advantage may be extended beyond A^ = 3. The number of selective 
rotations required for /g points to imperfections within selective pulses as one likely source of 
error. In particular, it must be noted that possible effects of scalar coupling evolution during 
application of selective rotations were ignored. Indeed, for the indirectly coupled realization 
of f/^° the total duration of all the selective rotations is comparable to 1/2 J21. To the best 
of our knowledge this issue has not been addressed satisfactorily. A further possible source 
of error are inhomogeneities in the RF magnetic fields used for selective rotations. 

To conclude, we have provided a comprehensive NMR realization of the A^ = 3 Deutsch- 
Jozsa algorithm, at the same time demonstrating quantum information processing via a 



linear configuration of couplings and indirect realizations of two-qubit gates. The use of 
appropriate SWAP gates allows for the extension of our method to quantum computation 
with larger numbers of qubits. 

This work was supported, in part, by the DARPA and the ONR. We would also like to 
thank Gary Sanders for useful discussion. 
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FIGURES 

FIG. 1. Deutsch-Jozsa algorithm. Initialization to the state: |a;) = {xn^i) ... \xo) where 
Xi £ {0, 1}. i?n(90) rotates each qubit through 90° about n = cos(j)x + sin(^y. For a single qubit 
Rn{0) '■= e~*"'°"^/^ where n.a = nx&x + %^2/ + nz^z. The /-controlled gate, Uf \x) := (— !)■' ^^^ \x), 
evaluates /. The expectation value of \x) {x\ on the output state answers the problem conclusively. 

FIG. 2. ^^C output spectra for alanine: (a) Fiducial spectrum, (b) /2 output spectrum, (c)/4 
output spectrum, (d) /y output spectrum, and (e) /g output spectrum. Insets provide enlargements 
of all antiphase and doubly antiphase multiplets. 



TABLES 

TABLE L Representatives of each class of admissible functions (N = 3) and the corresponding 
density operator after the function evaluation step, pf := UfR^y{90)pth- The form of the spin j 
multiplet depends on the term containing a factor of I^. Line intensities within a given multiplet 
are equal. I^ terms yield a multiplet with lines of equal phase. /*/^ yields an antiphase multiplet 
if Jij 7^ 0. If Jij = antiphase lines overlap and the entire multiplet disappears. I],lilz yields 
a doubly antiphase multiplet for spin j if Jjj, Jjk ^ 0. If Jij = or Jjk = the multiplet again 
disappears. Antiphase and doubly antiphase multiplets contain lines whose phases differ by vr. 
Multiplication by —1 corresponds to a phase shift of vr for each line in the corresponding multiplet. 

Representative, / pf 

Constant: 

J const = U ~J^x ~ ^x ~ ^x 

Balanced: 

/i = X2 +4 - 4 - 4 

h = X2®xi® xo +ll + il+i^ 

Ja = X2XX e Xo -2/^4^ - 2/1 /i + /° 

/s = X2X\ © X2 e Xo +2/J/i - 2/|/^ + /° 

/e = X2X1 © X2 © xi © Xo +2/x^i + 2/2 /^ + /° 

/y = X2X1 © XiXo © X2 © Xi +2/^^/] + 4/2/1/0 - 2/]/° 

/s = X2X1 © xixo © X2 +2/:^/i - 4/f 1^/° - 2/^/0 

/g = X2X1 © XiXo © X2X0 -4/2/1/0 - 4/2/1/0 - 4/2/1/0 
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Function: f1(x2x1x0) = x2 
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Deutsch - Jozsa Output 

Function: f2(x2x1x0) = x2 + x1 
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Deutsch - Jozsa Output 

Function: f3(x2x1x0) = x2 + x1 + xO 
Input: Thermal Equilibrium 
Readout: None 
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Deutsch - Jozsa Output 

Function: f4(x2x1x0) = x2x1 + xO 
Input: Thermal Equilibrium 
Readout: None 
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Deutsch - Jozsa Output 

Function: f5(x2x1x0) = x2x1 + x2 + xO 
Input: Thermal Equilibrium 
Readout: None 
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Deutsch - Jozsa Output 

Function: f6(x2x1x0) = x2x1 + x2 + x1 + xO 
Input: Thermal Equilibrium 
Readout: None 
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Deutsch - Jozsa Output 

Function: f7(x2x1x0) = x2x1 + x1xO + x2 + x1 
Input: Thermal Equilibrium 
Readout: None 
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Deutsch - Jozsa Output 

Function: f8(x2x1x0) = x2x1 + x1xO + x2 
Input: Thermal Equilibriunn 
Readout: None 
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Deutsch - Jozsa Output 

Function: f9(x2x1x0) = x2x1 + x2x0 + x1xO 
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Deutsch - Jozsa Output 

Function: f10(x2x1x0) = x2x1 + x2x0 + x1xO + x1 + xO 
Input: Thermal Equilibrium 
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